Experimental and First principle calculation of 
Co^Ni(i_^)Si solid solution structural stability 

J. Teyssier^, R. Viennois^, J. Salamin^, E. Giannini^ and 

D. van der MareP 

^Departement de Physique de la Matiere Condensee, Universite de Geneve, Quai 
Ernest- Ansermet 24, 1211 Geneve 4, Switzerland 



Abstract 

We report the investigation of the structural stabihty of Co(i_^)Ni2;Si monosih- 
cides for < X < 1. As CoSi crystalHzes in the FeSi-type structure (B20) and NiSi is 
stable in the MnP-type structure (B31), a complete set of samples has been synthe- 
sized and a systematic study of phase formation under different annealing conditions 
were carried out in order to understand the reason of such a structural transition 
when X goes from to 1. This study has revealed a limit in the solubility of Ni in 
CoSi B20 structure of about 17.5 at.% and of Co in NiSi B31 phase of about 13 at.%. 
For 0.35 < X < 0.74 both B20 and B31 phases are present in the sample at there 
respective limits of solubility. The temperature dependence of the magnetic suscep- 
tibility has also been measured revealing diamagnetic behaviors. Optimal structural 
parameters and phase stability of the solid solution have been investigated using 
self-consistent full-potential linearized augmented plane wave method (FP-LAPW) 
based on the density functional theory (DFT). This calculation well predicts the 
structural instability observed experimentally. 

Key words: A inter met allies, B crystal growth, C crystal structure, C phase 
transitions 



1 Introduction 

During the past years, the transition metal monosilicides MSi with B20 cubic 
structure have attracted lot of attentions due to their interesting and various 
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ground states. Notably, MnSi is an itinerant helimagnetic metal for T < 30K 
[1], FeSi is a paramagnetic Kondo insulator [2], CrSi is a Pauli paramagnetic 
metal [3,4] and CoSi is a diamagnetic metal [4,5,6]. Moreover, it should be 
noted that the solid solution Fea;Co(i_a^)Si also exhibits itinerant helimagnetic 
metallic behavior like MnSi for 0.4 < x < 0.9 (Tc = 60K for x = 0.6) although 
the two end-compounds FeSi and CoSi have no magnetic ordering [7,8,9,10]. 
On the other hand, NiSi is a diamagnetic metal which crystallizes in the 
B31 orthorhombic structure at room pressure [4,11]. Thus, it is interesting 
to study the solid solution Co(i_a^)Nia^Si for which no physical properties has 
been reported yet. Although the CoSi-NiSi phase diagram has already been 
investigated, very different limits of the solubility of Ni in the B 20 CoSi 
phase from 10 to 50% has been reported [12,13,14] depending on synthesis 
conditions and thermal treatments. On the other side of the phase diagram, 
only the Coo.2Nio.8Si in the B31 structure was synthesized[15]. 

The isothermal cross-section of the ternary phase diagram Co-Ni-Si at 800 
was reported by van Beek et al. [16]. Even though they investigated only the 
solid-solid equilibria at one given temperature they first reported the existence 
of a miscibility gap between CoSi and NiSi. 

In this paper, we report theoretical and experimental study of the structural, 
thermodynamic and magnetic properties of the whole range of Co(i_a;)Nia;Si 
compositions. The limits of solubility of Ni in B20 CoSi and Co in B31 NiSi are 
determined by combining X-ray diffraction, differential thermal analysis and 
quantitative EDX chemical analysis. We show that the stability ranges of the 
B31 and B20 structures are well predicted by total energy LDA calculations. 



2 Crystal structure 

CoSi crystallizes in the common FeSi-type B20 structure (P2i3 (198)) [17]. 
NiSi crystallizes in the MnP B31 structure (Pnma (62)) [18]. The existence of 
NiSi with the B20 structure has only been reported when part of the silicon 
is substituted with Al [19]. 




C) Metal ion Silicon ion 

Fig. 1. a) B20 structure of CoSi, b) B31 structure of NiSi 
The positions of both Metal and Si atoms in the B20 unit cell are (x,x,x). 
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and the 3 permutations of (x + — x^x). In CoSi, these values are found 
to be xco=0.14 and xsi=0.84:3 and the lattice parameter a=4.438 A [17]. 
In the B31 structure, the positions of both Metal and Si atoms are 

and For NiSi, these values are found to be 

xco=0. 00757, zco=0.18772,x^^=0. 32090 and z^^=0. 08168. Lattice parameters 
are a=5.194 A, b=3.323 A and c=5.629 A [18]. 

In the solid solutions of transition metal (Mn^Fe^Co) monosilicides, the 
cell parameter decreases, following the evolution of the ionic radius of the 
transition metal ion (Fig. 2). 
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Fig. 2. Evolution of the volume of the unit cell with chemical composition of 
monosilicides. Atomic radius values from Goldschmidt [20] are given as an indi- 
cation. 

A trend inversion, already reported by Watanabe et al. [13] (orange diamonds 
in Fig. 2), occurs between CoSi and NiSi for which the lattice parameter in- 
creases whereas the ionic radius of the metal ion is still decreasing. They no- 
ticed this evolution of the structure of Co(i_a;)Nia;Si solid solution until x ^ 0.2. 



3 Experimental 



All samples were synthesized using a home made arc furnace with a water 
cooled copper crucible, starting from 4N purity transition metals and 6N sili- 
con chunks. An annealing at 900^(7 from a minimum of 12 hours to a maximum 
of a week under high vacuum (about 5.10"'' mbars) is necessary to improve 
the crystalline order, increase the limits of solubility and decrease the devia- 
tion on X in solid solutions. This annealing temperature was chosen close to 
the melting point of NiSi (T = 979^C [21], T = 982^C [22]). Samples were 
processed every x = 0.1. 
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X-ray powder diffraction (XRD) was performed in a Philips PW1820 diffrac- 
tometer using the K^, radiation of a Cu tube (A = 1.5406 A). The XRD spectra 
were analyzed with a full pattern profile refinement method using the FuUprof 
program suite [23]. 

The stability range and melting temperature of each composition where mea- 
sured by DiflPerential Thermal Analysis (DTA) in a SETARAM TAG 24 ther- 
mal analyzer using AI2O3 as a reference. DTA measurements were performed 
under flowing Ar (0.6 1/h). About 50 mg of each sample were subjected twice 
to the same run at 5 ^C/min. The melting point of end compounds was iden- 
tified as the onset of the second heating endothermic peak. The Liquidus 
temperatures for solid solutions and mixed phased samples were assumed to 
be the offset of the second broad endothermic peak during the second heating 
run. The analysis of the chemical composition was carried out by Electron 
Dispersive X-ray spectroscopy (EDX) in a LEO 438VP electron microscope 
coupled to a Noran Pioneer X-ray detector, at a beam energy of 20 KeV. 
Quantification of elements was done on the K-lines using internal calibration. 



4 Calculation of structural stability 

The ah initio structural optimizations and total energy calculations were car- 
ried out using the density functional theory DFT method as implemented in 
the Quantum-ESPRESSO package [24]. The calculations were carried out us- 
ing an exchange-correlation functional by Perdew, Burke, and Ernzerhof [25] 
and Vanderbilt-type ultrasoft pseudopotentials [26]. Wave functions were ex- 
panded in a plane wave basis to a 80 Ry cutoff. A kinetic energy cutoff of 200 
Ry was used for the charge density. A 3x3x3 k-mesh has been used to sam- 
ple the Brillouin zone. First, the atomic positions and lattice constants were 
optimized for CoSi, Coo.5Nio.5Si and NiSi in both B20 and B31 crystal struc- 
ture. The starting point was the parameters extracted from X-ray refinements 
when available. For CoSi and Coo.5Nio.5Si in the B31 configuration we took 
the parameters of NiSi and for NiSi and Coo.5Nio.5Si in the B20 configuration, 
we took parameters of CoSi. 

In table 1, we report a summary of the lattice parameters obtained experi- 
mentally in this work from powder diffraction, from structural optimization 
and from the literature. 



The total energy of Co(i_^)Nia;Si was then calculated for x = 0, 0.25, 0.5, 0.75 
and 1. The calculation has been performed using the experimental parameters 
when available. When the cell parameters were not experimentally accessible. 
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CoSi opt. 


4.442 


0.143 
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0.3126 


0.0582 


CoSi [17] 
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0.843 
















NiSi exp. 








5.1818 


3.334 


5.619 


0.0069 


0.1876 


0.3122 


0.1041 


NiSi opt. 


4.515 


0.146 


0.846 


5.2972 


3.2549 


5.6997 


0.0068 


0.1884 


0.3181 


0.0786 


NiSi [18] 








5.194 


3.323 


5.629 


0.0076 


0.1877 


0.3209 


0.0817 


Coo.5Nio.5Si opt. 


4.475 


0.144 


0.846 


5.369 


3.029 


5.910 


0.0053 


0.1929 


0.3096 


0.0640 



Table 1 



Comparison of experimental and calculated lattice parameters and atomic positions 
for B20 and B31 structures for CoSi, Coo.5Nio.5Si and NiSi. exp. and opt. refer 
respectively to experimental and ab-initio optimized structures from this study. 
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we used the ones from structural optimization calculations. For the B20 struc- 
ture, the matching between calculated and experimental parameters is good 
in the range where both are available {0 < x < 0.35) (green open circles in 
Fig. 5b). This is not the case for the B31 structure where LDA overestimates 
the cell parameter (green closed circles in Fig. 5b). For the B31 structure, we 
used the composition dependence of the cell parameter obtained from calcula- 
tion scaled to match the experimental values in the range 0.74 < x < 1 (gray 
closed squares in Fig. 5b). 



The self consistent total energy calculation has been performed on a 10x10x10 
k-mesh of the irreducible Brillouin Zone. For x = 0.25, x = 0.5 and x = 0.75, 
we used a unit cell where respectively |, | and | Co atoms where substituted 
by Ni in both B20 and B31 structures. 

The realistic estimation of the value of the cell parameters allows us to calcu- 
late the total energy of all structures (real and hypothetic). In order to show 
how sensitive is the total energy to the deviation in cell parameters, we plotted 
in Fig. 3 the total energy difference on both "as-optimized" (open symbols) 
and "experimental" (closed symbols) structures. 




Fig. 3. (color online) Open and closed symbols correspond to the total energy calcu- 
lated from theoretical structural parameters and experimental refinements respec- 
tively. 

The total energy difference between B20 and B31 phases (Fig. 3) shows that 
the B20 structure is stable up to a concentration x 0.55. Above this con- 
centration, the B31 configuration saves energy. Even if this calculation cannot 
predict the limit of solubility, the domain of stability of the B20 (Ni in CoSi) 
is larger than the one of B31 phase (Co in NiSi). This agrees well with the ex- 
perimental observation of a much larger limit of miscibility of Ni in B20 CoSi 
than Co in B31 NiSi. In the following section, we discuss the experimental 
results. 
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5 Experimental results 



From systematic XRD, SEM and EDX analyses of the whole series of samples, 
we have identified three different composition ranges, characterized by three 
different sample morphologies when x is ranging from to 1: 

• a pure B20 solid solution; 

• a mixture of both B20 and B31 solid solutions; 

• a pure B31 solid solution. 

From XRD pattern refinement, the cell parameters and the ratio of the two 
different crystalline forms were extracted. The amount of the two different 
phases as a function of the nominal composition is plotted in Fig4. It clearly 
appears that below a certain level of substitution (Ni in B20-CoSi and Co in 
B31-NiSi) a good solubility is observed and samples are single phased. 



Fig. 4. B20 and B31 phase molar fractions present in the samples before and after 
annealing. The abscissa is the nominal composition. 

By combining the XRD with EDX analyses, it is possible to determine the 
exact composition and the cell parameters of each solid solution present in 
the multiphase mixture with 0.3 < x < 0.9. The amount of each phase as a 
function of real and nominal compositions, is shown in Fig. 5 a and Fig. 4, 
respectively. The evolution of the cell volume (i.e. average cell parameter) with 
the real composition determined from EDX is compared to LDA structural 
optimizations in Fig. 5b. 

In the limit of sensitivity of the XRD technique, the Co(i_a;)Nia;Si B20 phase 
is observed for x < Xi = 0.35 (17.5 at.%) and the B31 phase is observed for 
X > X2 = 0.74 (13 at.%). These two limits of solubility are displayed as vertical 
dashed lines in Fig. 5 and 4. We noticed a good agreement with van Beek et 
al. [16] who reported limits of solubility at 800^(7 of 12 at.% of Co in NiSi 
and 22 at.% of Ni in CoSi. We also learn from Fig. 4, in which these limits 
of solubility almost coincide with the limit where the samples present a single 
phase, that the solid solutions are stable up to the limit of solubility. 
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Fig. 5. a) Molar fraction of the B20 and B31 phases in the sample before and after 
annealing, b) Experimental and LDA predicted values of the cell volume for B20 
(closed symboles) and B31 (open symbols) structures. The vertical lines represent 
the two limites of solubility of Ni in CoSi and Co in NiSi. 



Within the stability range, LDA structural predictions for the B20 structure 
type agree perfectly with the experimental parameters. For the B31 struc- 
ture, energy minimum is found for a volume cell slightly higher than what is 
observed experimentally. A distorsion along the c axis is also observed. 

Between these two limits (for 0.35 < x < 0.74), the two solid solutions at the 
limit of solubility coexist. The Fig. 6 shows EDX Co and Ni maps in samples 
with nominal compositions x = 0.4, 0.5 and 0.6. 

The chemical composition and grains can be easily identified revealing sharp 
grain boundaries. The fact that cobalt rich grains are closed volumes comes 
from the fact that B20 cobalt rich structure has a higher melting point as we 
will see later. As a result, these grains nucleate and start growing first during 
cooling. B31 is growing at lower temperature in the free space between B20 
grains. 
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Fig. 6. The images correspond to Co (green) and Ni (red) EDX maps for different 
values of x in Co(i_3,)Nia;Si after annealing of the samples at 900^C during a week. 

This difference in melting temperatures was observed by thermal analysis for 
pure solid solution (Fig. 7b and d) and for multiphase samples (Fig. 7c). The 
evolution of the liquidus temperature with x is plotted in Fig. 7a. 
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Fig. 7. a) Evolution of the liquidus temperatures with x value in Co(i_3,)Nia;Si. 
Limits of solubility are displayed as dashed gray lines, b-d) thermal flow curves for 
respectively B20, mixed and B31 solid solutions. 



The melting temperature of NiSi T = 978^(7 agrees well with other recent 
DTA experiments (T = 979^C [21], T = 982^C [22]). Whereas no recent data 
are available for CoSi, the melting temperature T = 1446^(7 was found in the 
range of two old experimental reports (T = 1420^C [27] and T = 1460^C [28]). 
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These thermodynamic investigations first revealed that the liquidus temper- 
ature of these solid solutions decreases with x in the two region presenting 
a single phase. In the multiphase region (0.35 < x < 0.74), two samples 
have been measured (Coo.6Nio.4Si and Coo.4Nio.6Si) (Fig. 7c) pointing out a 
monovariant line at a temperature ^ 1040^(7 shown by points 6 and 8 in 
Fig. 7a. The second endothermic peak at higher temperature correspond to 
the liquidus. 

The magnetic susceptibility of some samples was measured from room tem- 
perature down to 4 K using a SQUID magnetometer showing diamagnetic 
behaviors over the whole temperature and composition range. 



6 Conclusions 



We have reported the study of the stability range of the two phases (B20 
and B31) forming in Co(i_a^)Nia;Si. The whole range of nominal composi- 
tion < X < 1 was investigated. Ab-initio calculation were used to opti- 
mize the structural parameters and predict the structural transition from to- 
tal energy computation. Systematic structural and thermodynamic studies of 
Co(i_a;)Nia;Si (0 < X < 1) revealed the existence of 3 different sample mor- 
phologies corresponding to 3 distinct regions of the CoSi-NiSi quasi-binary 
phase diagram: 

• A single phased B20 solid solution with a limit of solubility Xi = 0.35 

• A single phased B31 solid solution with a limit of solubility X2 = 0.74 

• and, between these two limit composition, a mixture with the two solid 
solutions at there limits of solubility. 

In the mixed phase range, a monovariant line at a temperature of 1040 
was found. 
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